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Abstract The GOES X28 flare of 4 November 2003 was the largest ever recorded in its 
class. It produced the first evidence for two spectrally separated emission components, one 
at microwaves and the other in the THz range of frequencies. We analyzed the pre-flare phase 
of this large flare, twenty minutes before the onset of the major impulsive burst. This period 
is characterized by unusual activity in X-rays, sub-THz frequencies, Ha, and microwaves. 
The CME onset occurred before the onset of the large burst by about 6 min. It was preceded 
by pulsations of 3 - 5 s periods at sub-THz frequencies together with X-ray and microwave 
enhancements. The sub-THz pulsations faded out as impulsive bursts were detected at 100- 
300 keV and 7 GHz, close to the time of the first Ha brightening and the CME onset. The 
activities detected prior to and at the CME onset were located nearly 2 arcmin south of the 
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following large flare, suggesting they were separate events. This unusual activity brings new 
clues to understanding the complex energy buildup mechanisms prior to the CME onset, 
occurring at a distinct location and well before the major flare that exploded afterwards. 

Keywords Coronal mass ejections (CMEs) ■ CME precursors • Solar flares • Sub-THz 
emissions • X-ray emissions 


1. Introduction 

The relationships between coronal mass ejections (CMEs) and flares remain unclear, mainly 
because there are a variety of situations that prevent a single explanation (see the review 
by Mittal and Narain (2010) and references therein). The initiation phase of CME events 
is not well known in most cases, because most coronagraphs cannot observe the source 
region in the solar atmosphere at lower altitudes. However, the analysis of a few events 
at lower altitudes in the solar atmosphere has suggested a two-phase CME launch process 
(Zhang et al, 2001). The initiation phase starts at slow speeds (tens of kms -1 ), tens of 
minutes before the acceleration (or onset) phase, exhibiting higher speeds (hundreds to a 
few thousand kms -1 ). The onset phase begins at altitudes lower than those in available 
coronagraph observations. 

There are indications that the CME onset phase is associated with weak soft X-ray en- 
hancements, sometimes preceding larger X-ray bursts (Datlowe, Elcan, and Hudson, 1974; 
Harrison et al, 1985; Gopalswamy et al, 2006; Schwenn et al, 2006). For a number of 
events analyzed it has been found that there were sub-THz pulsations associated with the 
CME onset phase (Kaufmann et al, 2003). On the other hand, there are CMEs occurring 
without any obvious flare connection (Webb, 2000). There have been studies suggesting that 
solar flares are a consequence of CMEs rather than causing them. Namely, the expansion of 
a CME loop may subsequently trigger magnetic reconnection, giving rise to a flare after- 
wards. This view is consistent with a general tendency that the CME onset precedes the 
flare (Harrison et al, 1985; Harrison, 1986; Kahler, 1992). On the other hand it has been 
suggested that both the CME initiation and the pre-flare phases involve ejection and heat- 
ing of material by magnetic reconnection, and therefore they may all be part of a common 
dynamic instability process (Ohyama and Shibata 1997, 1998). 

On 4 November 2003, the X-ray intensity of the flare started to increase at about 
19:32 UT, and peaked between 19:40- 19:50 UT. The major microwave and sub-THz inten- 
sities rose at about 19:40 UT (Kaufmann et al, 2004). These start times were derived as the 
time when intensities rose and became detectable, and therefore, they contain an uncertainty 
of several minutes. This flare is known as the largest soft X-ray burst ever recorded (extrap- 
olated GOES class X28 to X48, after Woods et al. (2004), Thomson, Rodger, and Dowden 
(2004), and Brodrick, Tingay, and Wieringa (2005)) and for producing the first evidence 
for a THz emission component spectrally distinct from the well-known microwave emission 
(Kaufmann et al, 2004). The Large Angle and Spectrometric Coronagraph (LASCO) in- 
strument (Brueckner et al, 1995) onboard the Solar and Heliospheric Observatory (SOHO) 
detected a large coronal mass ejection (CME) associated with this burst. We find that the 
CME onset phase, derived from the limb crossing time after a linear extrapolation back from 
the height-time coronagraph observations to the solar surface, occurred considerably before 
the onset time of this large burst. This extrapolation is shown in Figure 1 , suggesting that 
the CME onset phase began at about 19:34 UT, with ail min uncertainty (assuming the 
validity of the linear extrapolation). The figure also shows three data points from the NASA 
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Figure 1 Radial distances of CME features above the solar limb, derived from LASCO C2 (filled triangles) 
and C3 (filled circles) instruments. The limb crossing time by a linear extrapolation suggests an onset time at 
about 19:34 UT on 4 November 2003, close to the first Ha brightening observed by BBSO (filled diamonds, 
see Figure 3). The data from NASA CME catalogue (http://cdaw.gsfc.nasa.gov/CME_list/) have been added, 
with two positions obtained from C3 (open diamonds) and one from C2 (open square), exhibiting a similar 
onset time, within the uncertainty of a few minutes, assuming the validity of the linear extrapolation. 
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Figure 2 The first Ha brightening observed by Big Bear Solar Observatory around 19:32 UT on 4 November 
2003. The crosses on the right-most panel indicate the source positions of the sub-THz pulsation at about 
19:30 UT and the great flare occurring after 19:40 UT. 


CME catalogue (http://cdaw.gsfc.nasa.gov/CME_list/) (which gives a limb crossing time of 
19:42 UT or 19:36 UT for rough linear or quadratic extrapolations, respectively). The onset 
time defined by these extrapolations, however, has an uncertainty due to the limited field of 
view of the LASCO C2 coronagraph at lower altitudes (< 10 5 km). As noted earlier, it is 
possible that at altitudes lower than 10 5 km an initiation phase, rising with much smaller 
speed, might have started tens of minutes earlier (Zhang et al., 2001). The presence of the 
initiation phase could not be directly observed in the LASCO C2 and C3 data analyzed here. 
On the other hand, the Ha movie obtained by Big Bear Solar Observatory (BBSO, see http:// 
www.bbso.njit.edu/) shows the very first brightening at about 19:33 - 19:35 UT, close to the 
suggested CME onset time (see the four successive Ha images in Ligure 2). 

The Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) space exper- 
iment (Lin et al, 2002) obtained high quality X-ray observations just before the large 
burst of 4 November 2003. These data were investigated together with observations ob- 
tained at sub-THz frequencies by the Solar Submillimeter Telescope (SST; Kaufmann et 
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al., 2008), in microwaves obtained by the solar array at Owens Valley Radio Observatory 
(OVRO; Gary and Hurford, 1999), and the Itapetinga 7 GHz polarimeter (Kaufmann, 1971; 
Correia, Kaufmann, and Melnikov, 1999). Our investigation reveals an unusual pulsating 
and bursting activity, together with a slow, varying enhancement in X-rays, preceding the 
CME onset. We discuss a possible association between the early CME initiation phase and 
the major burst. 


2. Multi-wavelength Observations 

In Figure 3 we show the time profiles describing the solar activity 20 min prior to the large 
burst. The dashed vertical line indicates the approximate CME onset time. GOES 1 - 8 A and 
RHESSI X-ray time profiles are displayed at the top. The GOES soft X-ray flux appeared 
to rise after about 19:29 UT, nearly the same time as the rise in RHESSI 6-10 keV X- 
rays. The RHESSI 25-50 keV X-ray flux began to rapidly increase after 19:32 UT. This 
rapid increase is preceded by two slow oscillations at about 19:27:30 UT and 19:30 UT. 
The sudden intensity drop just before 19:34 UT (19:33:48 UT) in the 6- 10 and 25-50 keV 
bands is an artifact caused by the automatic insertion of X-ray attenuators. The enhanced 
count rate in these bands from 19:38:00 to 19:38:28 UT was caused by a brief removal of 
the attenuators. The RHESSI 100-300 keV flux exhibits one impulsive burst peaking at 
about 19:34 UT, which is nearly coincident with a microwave impulsive burst (7 GHz flux 
shown in the bottom of the figure). Both are nearly coincident with the derived CME onset 
time. These events were followed by another small impulsive burst, detected in the two 
higher energy RHESSI X-ray bands only, at about 19:36.3 UT. RHESSI was in spacecraft 
nighttime after 19:39 UT. 

The bottom panel of Figure 3 presents the SST data of beams 2 (0.2 THz) and 5 (0.4 THz) 
in units of scintillation index (SI). This index is used to detect the occurrence of pulsations, 
to facilitate their identification and extraction from the data. It is a normalized, dimensionless 
magnitude defined as the standard deviation (sigma) of fluctuating data sampled every 40 ms 
relative to the signal average taken every 3 s, divided by the difference of the signals obtained 
by the antennas pointing to the Sun and pointing to the sky. The SI data displayed in Figure 3 
have been smoothed by a running mean over 10 s. The pulsations were enhanced around 
19:25 UT, with SI exhibiting slow modulations of the order of 1.8 min at both 0.2 and 
0.4 THz. The pulsations clearly faded away after the occurrence of the hard X-ray and 
microwave bursts at about the time of the CME onset or at the beginning of an acceleration 
phase around 19:34 UT. 

The time profile of 7 GHz total flux (R + L) is also shown in the bottom panel of Figure 3. 
A small precursor started at about 19:28 UT, together with an increase in the RHESSI X- 
rays. Two 25-50 keV X-ray enhancements at 19:27:30 and 19:30 UT were observed as 
the sub-THz pulsations became more pronounced. These were followed by the microwave 
impulsive burst, and almost simultaneously, by the high-energy RHESSI impulsive event. 
Both events occurred close to the derived CME onset time (19:35 UT). The microwave 
spectrum of the impulsive burst had a maximum at about 8 GHz (Figure 4). No clear sub- 
THz counterpart was observed by SST. 

Figure 5 shows the sub-THz pulsations during a portion of the period of enhanced scintil- 
lation index (Figure 3), in the 30 s interval 19:30:30- 19:31:00 UT. The two plots represent 
SST data of 40 ms sampling from beams 2 (0.2 THz) and 5 (0.4 THz), without any running 
mean smoothing. Despite the much noisier 0.4 THz data, there is a nearly one-to-one corre- 
spondence of quasi-periodic time structures, displaying about 14 time structures in 30 s. 
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Figure 3 Multi-wavelength time 
profiles observed on 4 November 
2003 before the onset of the large 
burst at 19:40 UT. As the GOES 
and RHESSI soft X-rays 
increased (upper two plots in the 
top panel), there was an 
enhancement of the scintillation 
indices at two sub-THz 
frequencies (bottom panel). 
Impulsive bursts were observed 
in hard X-rays and microwaves at 
about the CME onset time. The 
sudden drops in RHESSI 6-10 
and 25 - 50 keV signals at about 
19:34 UT are instrumental 
effects. 
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Figure 4 Spectrum of the 
microwave impulsive burst that 
occurred at about 19:34 UT, 
obtained by the Owens Valley 
Solar Array. 
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Figure 5 Enlarged SST flux variations showing pulsations observed at 19:30:30- 19:31:00 UT. There is a 
nearly one-to-one correspondence between the pulses at 0.2 (top panel) and 0.4 (bottom panel) THz, for about 
14 time structures. 


3. Spatial Positions 

The SST has an arrangement of four 0.2 THz and two 0.4 THz feed horns at the focal 
plane of a 1.5 m diameter antenna (Kaufmann et al., 2008). The corresponding half-power 
beamwidths are 4 arcmin and 2 arcmin, at 0.2 and 0.4 THz, respectively. Three of the 
0.2 THz beams (2, 3, and 4) are arranged in a triangular configuration and partially overlap 
with each other with their axes separated by the half-power beamwidth. One of the 0.4 THz 
beams (beam 5) points at the center of the triangle. This cluster of four beams tracks the 
solar active center. The other two beams 1 (0.2 THz) and 6 (0.4 THz) are co-aligned and 
displaced by 8 arcmin from the cluster of the four beams (i.e., beams 2, 3, and 4 at 0.2 THz 
and beam 5 at 0.4 THz). A comparison of partially overlapping beams at 0.2 THz allows the 
determination of the burst position in space (Gimenez de Castro et al., 1999). 

In Figure 6 we show the SST source positions with the crosses on a 195 A EUV image 
from SOHO EIT. The size of the crosses indicates the approximate uncertainty of the po- 
sition determinations which is about ± 25 arcsec. The source position of the pulsations at 
19:30 UT is about 2 arcmin south of the position of the major burst detected nearly 10 min 
later (19:40- 19:50 UT) and very close to the intensity peaks in RHESSI images at 12- 
20 keV and 25 - 50 keV, as shown in Figure 6. At an earlier time (19:29:12 UT) the RHESSI 
source was further north and between the two SST sources indicated by crosses in Figure 6. 
The RHESSI source location matched closely the source position of the SST pulsations as 
they became more pronounced, at about the time of the impulsive bursts in microwaves and 
X-rays, and the CME onset time (i.e., at 19:33:50 UT). These positions are also close to the 
first chromospheric Ha brightening shown in Figure 2. 


4. Discussion 

The following remarkable features were found for this event, in the pre-flare phase of the 
large X28 flare: 
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Figure 6 The RHESSI 12-20 and 25-50 keV X-ray images and SST source positions overlaid on an EIT 
195 A solar disk image. Crosses show the SST source positions of the pulsations at 19:30 UT (red), and the 
large burst that occurred later at 19:44 UT (cyan), with uncertainties of ±25 arcsec. The RHESSI contours 
are shown in green (12-20 keV) and blue (25-50 keV). Panel (a), obtained earlier at 19:29:12 UT, shows 
the X-ray sources located between the two SST crosses. Panel (b), obtained at 19:33:50 UT during the initial 
impulsive burst and at the derived CME onset time, indicates that the X-ray sources coincide in position with 
the SST pulsating source, nearly 2 arcmin south of the position of the large burst that began around 19:40 UT. 


i) The pulsations of emission at sub-THz frequencies started at least 10 min before the 
deduced CME onset time. The pulsations exhibited short periods of 3 - 5 s, and were 
modulated by a slower component with periods of about 1-2 min. 

ii) The sub-THz pulsations became more pronounced and the RHESSI 25 - 50 keV X-rays 
showed superimposed time structures as the intensity rose in RHESSI 6-10 keV, GOES 
soft X-rays, and in microwaves. 

iii) As the sub-THz pulsations faded out an impulsive microwave burst occurred together 
with a RHESSI hard X-ray impulsive burst, at about the CME onset time and the time 
of the first Ha brightening. 

iv) The sub-THz, RHESSI X-ray, and BBSO Ha activities occurred before and at about the 
CME onset time at a different location from that at the large burst which occurred nearly 
six minutes after the CME onset time. 

This detailed description reveals complex energy build up activity in X-rays, Ha, sub- 
THz wavelengths, and microwaves. The sub-THz pulsations became more pronounced after 
about 19:24 UT, together with the first soft X-ray rise observed by GOES and RHESSI. They 
probably represent physical processes corresponding to the CME initiation phase (Zhang et 
al., 2001) and are related to the subsequent CME onset phenomena, but apparently without 
a clear connection to the large flare that occurred afterwards. The whole picture fits well 
the evidence that the CME onset most commonly precedes the flare (Harrison, 1986; Lin, 
2005). 

The emission mechanism responsible for the sub-THz radiation has not been estab- 
lished. Two possibilities are free-free radiation from hot, thermal plasma and gyrosyn- 
chrotron from high-energy, non-thermal electrons in a strong magnetic field (Zirin, 1988; 
Zirin, MacKinnon, and Mackenna-Lawlor, 1991; and references therein). Either mechanism, 
however, must be optically thick, since the intensity at 0.4 THz is greater than at 0.2 THz. 
Below we consider implications of each mechanism. 
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The X-ray events prior to the CME onset time seem to be due to thermal emis- 
sion, favoring the suggested mechanism of reconnection with plasma heating and a pos- 
sible plasmoid ejection leading to the main CME acceleration phase (Ohyama and Shi- 
bata 1997, 1998). The sub-THz pulsations, with time scales of a few seconds, might be 
a response of the magnetically confined and heated thermal plasma to magnetohydrody- 
namic oscillations associated with resonant Alfven waves (for example Nitta et al., 1991; 
Fleishman, Bastian, and Gary, 2008; Kislyakova et al., 2011). These modulations have an 
important role in producing the RHESSI 25 - 50 keV time structures, in the time interval 
of the first soft X-ray enhancement, which represent a response to reconnection episodes 
at different heights, with possible plasmoid ejections adding to the overall CME evolution 
in space (Ohyama and Shibata 1998, 1998; Nishizuka et al., 2010). The fading out of the 
pulsations at the CME onset time might be related to the ultimate release of heated plasma, 
leaving only a relatively small amount of material, not sufficient to respond to the magnetic 
field modulations. 

The processes occurring during the CME initiation phase and prior to the onset time 
seem to include the clustering of rapid and successive pulses at sub-THz frequencies. This 
can be reconciled with other non-thermal mechanisms closer to or at the solar surface, such 
as the superposition of multiple or fragmented energy releases (Biesecker and Thompson, 
2000). This process produces gyrosynchrotron emissions, as indicated by a number of flare 
diagnostics and model descriptions (Frost, 1969; van Beek, de Feiter, and de Jager, 1974; 
Kaufmann et al., 1980; Kaufmann, 1985; Lu and Hamilton, 1991; Zirker and Cleveland, 
1993; MacKinnon, MacPherson, and Vlahos, 1996). As the magnetic complexity grows in 
an active region, conditions are reached before a major large-scale instability sets in, induc- 
ing multiple and faster instabilities on smaller scale magnetic arcades or fluxules (Sturrock 
and Uchida, 1981; Sturrock et al., 1984; Sturrock, 1986; and references therein). The en- 
ergy released by such elementary instabilities produces nanoflares or, when clustered in 
a larger number in a short time, brighter flares. Models based on larger spatial scale struc- 
tures describe CMEs arising from the disruption of certain coronal streamer configurations 
by flux -rope instabilities at the photospheric level (Wu et al . , 2000), which can be reconciled 
with similar scenarios with multiple and complex magnetic structures. 

Quasi-periodic pulsations (QPP) with time scales of seconds to tens of seconds obtained 
at short microwave wavelengths (17 and 34 GHz) with the Nobeyama Radioheliograph have 
been studied by Inglis, Nakariakov, and Melnikov (2008). Two groups of generating mech- 
anisms have been proposed (Melnikov et al., 2005; Nakariakov and Melnikov, 2009). The 
former assumes repetitive regions of flaring energy releases by magnetic reconnection or 
other means, and the latter assumes coherent MHD oscillations of all regions contributing 
to the flare emission. We may extend this interpretation to the sub-THz pulsations. The sam- 
ple shown in Figure 5 indicates that the pulses are in phase at the two frequencies with a 
reduction in amplitude for the higher 0.4 THz frequency. This may suggest optically thick 
emission from very dense plasma (> 10 12 cm -3 ). Adopting the global sausage mode of oscil- 
lations in a magnetic tube (Melnikov etal., 2005), the sub-THz radiation possibly originates 
from the inner, highly dense sections of the tube. 

Finally, the assumption of multiple energy releases as the origin of the pulsations can 
also be reconciled with recent results obtained from simulations of far infrared and sub-THz 
thermal continuum emission as a response of denser plasma in the solar atmosphere to time- 
variable heating by electron beams (Kasparova et al, 2009). They show that the continuum 
intensities can vary on the same time scale as short-duration variations in the electron beam 
flux, from sub-second to longer durations. 
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5. Conclusions 

We analyzed in detail the activities 20 min before the large burst of the 4 November 2003 
flare with instruments at multiple wavelengths: X-rays by RHESSI, LASCO coronagraphs 
on the SOHO spacecraft, Ha by BBSO, two sub-THz frequencies by SST, and high sen- 
sitivity Itapetinga polarimeter at 7 GHz, and at multiple microwave frequencies by OVRO. 
The first enhancements at soft X-rays and microwaves were observed nearly 10 min before 
the CME onset. We also found the sub-THz pulsations nearly 6 min before the onset of the 
largest X-class flare. 

The CME onset time coincided with an impulsive burst in hard X-rays and microwaves 
with moderate intensity, nearly 6 min before the large burst onset. The CME onset time 
was preceded by pulsations of 3 - 5 s periods at sub-THz frequencies that began nearly 
10 min earlier, together with soft X-ray and 7 GHz enhancements and superimposed 25- 
50 keV X-ray time structures. These unusual phenomena preceding the CME onset strongly 
indicate that the CME material started to be lifted up at least 10 min before a sudden energy 
release occurs and the CME is accelerated strongly upward (Ohyama and Shibata 1998, 
1998; Nishizuka et al., 2010). 

The position of the X-ray source corresponded to the approximate position of the sub- 
THz pulsations before and at the CME onset time, when the first Ha brightening appeared, 
all displaced by about 2 arcmin south from the position of the large flare that occurred 
later, suggesting they were distinct events. The sub-THz pulsations faded out as impulsive 
bursts were observed in X-rays and microwaves at about the same time as the CME onset. 
The complex, unusual activity before and at the CME onset time might be representative of 
the initiation phase preceding the sudden acceleration of the ejected material (Zhang et al., 
2001). 

The physical mechanism producing the sub-THz pulsations is still unclear. To improve 
our understanding, more observational and theoretical investigations are required. Exist- 
ing developments are potentially helpful in the search for explanations. For example, we 
may extend to sub-THz frequencies the adoption of two groups of mechanisms proposed 
to explain short microwave (i.e., 17 and 34 GHz) QPP pulsations (Melnikov et al, 2005; 
Nakariakov and Melnikov, 2009). One assumes multiple flaring energy releases, and an- 
other assumes MHD-modulated plasma emission. The response to MHD sausage-mode 
modulation of a confined plasma also has an important role in producing the RHESSI 
25 - 50 keV time structures. This in turn represents a response to reconnection episodes 
at different heights, with plasmoid ejections adding to the overall CME evolution in space 
(Ohyama and Shibata 1998, 1998; Nishizuka et al, 2010). The sub-THz pulsations dur- 
ing the CME initiation phase can also be reconciled with suggested mechanisms closer to 
or at the solar surface, such as the superposition of multiple or fragmented energy release 
events. As the magnetic complexity grows in an active region, even before a major large- 
scale instability sets in, multiple and faster instabilities may appear in smaller scale mag- 
netic arcades or fluxules (for example, Sturrock and Uchida, 1981; Sturrock et al, 1984; 
Sturrock, 1986). Recent results from simulations show rapid pulsations in far infrared and 
sub-THz thermal continuum emission as the response of denser plasma in the solar atmo- 
sphere to time-variable heating by electron beams with sub-second or longer durations (Kas- 
parova et al., 2009). 

A number of known CME-related solar activities were confirmed in this event, namely 
the onset of soft X-rays, the pronounced sub-THz pulsations before the CME onset time, and 
the close time coincidence between the CME onset and the earlier, smaller impulsive burst. 
New features were found here, such as the slow modulation (tens of seconds) of sub-THz 
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pulsations (seconds), which are somewhat related to enhancements in soft X-rays and time 
structures in harder RHESSI X-rays. The pulsations seemed to fade away once the CME was 
launched. The relative positions in space of the X-ray source, the sub-THz pulsations, and 
first Ha brightening before and at the CME onset, displaced by about 2 arcmin south from 
the position of the large flare that occurred later, suggests that they were distinct events. 

The overall scenario with respect to the large burst that occurred later can be well de- 
scribed by the qualitative picture (Harrison, 1986) of a CME loop expansion triggering a 
major flare away from the centerline of material ejection, tens of minutes after its launch 
time. 

New patterns associated with the CME occurrence were obtained from the combined 
analysis of the low-level X-ray activity observed by RHESSI together with sub-THz and 
microwave pulsations. Analysis of a larger number of events, for which there is a similar 
combination of data available, should reveal new clues about their physical nature and the 
links between them, which are not currently well understood. 
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